Regularizing the Divergent Structure of the Light-Front Currents 



Bernard L. G. Bakker^, Ho-Meoyng Choi^'^, and Chueng-Ryong Ji^ 
^ Department of Physics and Astrophysics, Vrije Universiteit, De Boelelaan 1081, NL-1081 HV 

Amsterdam, The Netherlands 
^ Department of Physics, North Carolina State University, Raleigh, NC 27695-8202, USA 
^ Physics Department, Carnegie Mellon University, Pittsburgh, PA 15213, USA 



Abstract 



The divergences appearing in the 3+1 dimensional fermion-loop calcu- 
lations are often regulated by smearing the vertices in a covariant manner. 
Performing a parallel light-front calculation, we corroborate the similarity be- 
tween the vertex-smearing-technique and the Pauli-Villars regularization. In 
the light-front calculation of the electromagnetic meson current, we find that 
the persistent end-point singularity that appears in the case of point vertices 
is removed even if the smeared- vertex is taken to the limit of the point-vertex. 
Recapitulating the current conservation, we substantiate the finiteness of both 
valence and non-valence contributions in all components of the current with 
the regularized bound-state vertex. However, we stress that each contribu- 
tion, valence or non-valence, depends on the reference-frame even though the 
sum is always frame-independent. The numerical taxonomy of each contribu- 
tion including the instantaneous contribution and the zero-mode contribution 
is presented in the n, K, and D-meson form factors. 
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I. INTRODUCTION 

With the recent advances in the Hamiltonian renormahzation program, Light-Front Dy- 
namics (LFD) appears to be a promising technique to impose the relativistic treatment of 
hadrons. In LFD a Fock-space expansion of bound states is made to handle the relativis- 
tic many-body effects in a consistent way |T|. The wave function ipn{xi, ki±, Xi) describes 
the component with n constituents, with longitudinal momentum fraction Xi, perpendicu- 
lar momentum ki± and helicity Aj, z = 1, . . . , n. It is the aim of LFD to determine those 
wave functions and use them in conjunction with hard scattering amplitudes to describe the 
properties of hadrons and their response to electroweak probes. Important steps were taken 
towards a realization of this goal 0. However, at present there are no realistic results avail- 
able for wave functions of hadrons based on QCD alone. In order to calculate the response 
of hadrons to external probes, one might resort to the use of model wave functions. The 
variational principle enabled the solution of a QCD-motivated effective Hamiltonian and the 
constructed LF quark-model provided a good description of the available experimental data 
spanning various meson properties 0. The same reasons that make LFD so attractive to 
solve bound-state problems in field theory make it also useful for a relativistic description 
of nuclear systems. Presently, it is realized that a parametrization of nuclear reactions in 
terms of non-relativistic wave functions must fail. LF methods have the advantage that 
they are formally similar to time-ordered many-body theories, yet provide relativistically 
invariant observables. Furthermore, as far as the concerned amplitude is unconditionally (or 
absolutely) convergent, the LF Hamiltonian approach must yield the same result as the one 
obtained by a covariant Feynman approach. 

However, not all is well. As we have recently shown the amplitudes that are only 
conditionally convergent must be treated with care. A case in point discussed in our work @] 
was the calculation of a current matrix element in quantum field theory. A typical amplitude 
is given by the triangle diagram. One encounters this diagram e.g. when computing the 
meson form factor (see Fig. |I](a)). The vertices denoted by $ are coupling constants in 
covariant perturbation theory. The hard scattering process is the absorption of a photon of 
momentum g by a (anti-)quark. In the LFD approach the covariant amplitude is replaced by 
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a series of LF time-ordered diagrams. In the case of the triangle diagram they are depicted in 
Figs. |l](b) and 1(c). The first (Fig. 0(b)) of these two diagrams is easily interpreted in terms 
of the LF wave functions \Ef. However, the other diagram (Fig. |l|(c)) has one vertex that can 
again be written in the same way as before, but it contains also another vertex, denoted by 

that cannot be written as a LF wave function. The necessity of this new element in 
LFD has also been discussed in the cases of semileptonic meson decays and deeply virtual 
Compton scattering . One may call and the vertices of leading order and non- leading 
order, respectively, in the sense that the leading order vertex corresponds to the lowest Fock- 
state whereas the non-leading order vertex takes into account the higher Fock-states. The 
diagrams with \E'' are designated as nonvalence diagrams while those with vertices of type 
\E' only are designated as valence diagrams. In order to obtain the invariant form factor, 
the two LF form factors must be added. It depends on the situation whether one can limit 
oneself to a single component of the current J'^ to extract the invariant objects. For the 
electromagnetic current of a spin-0 particle any single component would suffice to extract 
the unique form factor. On the other hand, in situations like semileptonic pseudoscalar 
meson decay, which involves two independent form factors, or the electromagnetic current 
of a (axial-)vector particle that is described by three independent form factors, one must 
use information from several current components to determine the invariant amplitudes. 

Earlier, we presented an analysis of contributions from the nonvalence diagrams 0. 
We constructed both leading and non-leading order vertices using pointlike covariant ones. 
The model that we used was essentially an extension of Mankiewicz and Sawicki's (1 + 1)- 
dimensional quantum field theory model [0 , which was later reinvestigated by several others 
p|^l2|. While their model |^ was a simple 1+1 dim. scalar field theory, it included a 



binding effect of the two-body bound-state. Indeed, in Ref. g], the relativistic two-body 
bound-state form factor was discussed in the full range of the binding energy. The starting 
model wave function was the solution of the covariant Bethe-Salpeter equation in the ladder 
approximation with a relativistic version of the contact interaction ||^ . The covariant model 
wave function was a product of two free single particle propagators, the overall momentum- 
conserving Dirac delta function, and a constant vertex function. Consequently, all our 
form factor calculations were various ways of evaluating the Feynman triangle diagram in 



quantum field theory. As pointed out in Ref. |13|, however, the elastic electromagnetic form 



factors of a bound-state computed from the triangle diagram and from the Hamiltonian 
front-form dynamics are the same |T^. Since our aim was to analyze the taxonomy of 



the triangle diagram, we didn't choose any particular gauge for the electromagnetic gauge 
field but presented the equivalence of the physical form factor F(g^) in any choice of the 
electromagnetic gauge Q or in any choice of the current component. 

Our conditionally convergent example of LFD with a fermion-loop showed that the bad 
component of the current, J~ , with spin-1/2 constituents exhibits a persistent end-point 
singularity in the contribution from the nonvalence diagram ||^. However, the calculation 
carried out so far was semi-realistic as the model was 1+1 dimensional and only a point- 
vertex was considered. 



^For a recent advocacy of the anti-light-cone gauge, A =0, see M. Morara, R. Soldati, and G. 
McCartor, AIP Conference Proceedings 494, 284-290 (1999). 
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In the present work, we extend our analysis of the fermion-loop to the case of 3+1 
dimensions. In 3+1 dimensions both the covariant and the LF calculations are divergent 
and the model without any smeared vertex for the fermion loop is not well defined. In the 



recent literature ||T5|,|16[, the fermion-loop was regulated by smearing the qq bound-state 
vertex in a covariant manner. However, the vertex function was not symmetric in the four 
momenta of the constituent quarks and could hardly be considered a realistic approximation 
of a qq bound state. It was regarded only as a convenient cutoff prescription which makes the 
one-loop integrals finite [^,0. Furthermore, the calculation of the meson decay constant 
reveals that the end-point singularity is not completely cancelled by such an asymmetric 
choice of the vertex function. (See the next section for more details.) We show in this work 
that the fermion-loop can also be regulated by taking a non-local gauge-boson vertex. With 
this method satisfying the Ward-Takahashi (W-T) identity we found the complete 



cancellation of the end-point singularity not only in the electromagnetic form factor but 
also in the decay constant. The non-local gauge-boson vertex remedies also the conceptual 
difficulty associated with the asymmetric way of treating q and q in the previous calculations 

EH- 

Nevertheless, one should distinguish the bound-state from the confined-state. In this 
work, we are treating the mesons as bound-states rather than confined-states, because we 
do not yet know how to make a covariant regularization for the confined-state. Thus, our 
emphasis here is the inclusion of the non-leading order vertex rather than the model-building 
of a realistic meson wavefunction. 

We have performed the LF calculation in parallel to the covariant Feynman calculation. 
Our light-front results entail the similarity between the vertex-smearing-technique (either 
for the bound-state vertex or the gauge-boson vertex) and the Pauli-Villars regularization. 
For the bosonic loop calculation, the two methods turn out to be identical. However, for 
the fermionic loop calculation, the vertex-smearing-technique shares only the same structure 
of the denominators with the Pauli-Villars regularization. Using the gauge-boson vertex- 
smearing technique, we found that the persistent end-point singularity is removed even if 
the smeared-vertex is taken to the limit of the point-vertex. 

A significant entity of our work is the taxonomy of valence and non- valence contributions 
substantiating the finiteness of each contribution when the gauge-boson vertex is regulated. 
Our results satisfy current conservation. However, we note that each contribution individu- 
ally depends on the reference-frame even though the sum is always frame- independent. We 
thus elaborate the frame- dependence of individual contributions. Also, the zero-mode con- 
tribution should be distinguished from the instantaneous contribution. For the numerical 
estimates of physical observables, we present the electromagnetic form factors of the vr, K, 
and D mesons. 

In the next Section (Section |l|), we present both the covariant Feynman calculations and 
the LF calculations using the LF energy integration for the electromagnetic form factors of a 
pseudoscalar meson with spin-1/2 constituents. Section |T| contains the numerical taxonomy 
of both the valence and nonvalence diagrams to the electromagnetic form factors of tt, K 
and D-mesons. The conclusion and discussion follow in Section The general formula 
including unequal mass cases are summarized in Appendix A and the analytic behavior of 
the valence and nonvalence contributions at ^ limit in various frames is summarized 
in Appendix B. 
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II. CALCULATIONS 



The electromagnetic form factors can be extracted from the matrix elements of the 
current 

{p'\j^\p) = zem{p"+pnn<i'), (2.1) 

where is the charge of the meson and = {p' — pY is the square of the four momentum 
transfer. If one uses the plus-component, = + J^, the LF calculation gives two finite 
contributions, the LF valence form factor and the LF nonvalence form factor, that add up 
to the covariant result, as expected. The importance of the nonvalence contribution varies 
strongly with the momentum transfer and depends sensitively on the binding energy of the 
meson. For small values of and small binding energy, the valence part is dominant, but 
elsewhere the nonvalence diagram is essential for agreement between the LF calculation and 
the covariant results. 

The form factor can also be extracted from the minus-component of the current, 
J~ = — J^. Covariance guarantees that it makes no difference whether the form fac- 
tor is determined using the plus or the minus current matrix element. As LFD is not 
manifestly covariant, it may happen that J~ leads to a form factor different from the one 
determined using J"*". As we showed in Ref. 0], the matrix element of J~ diverges even in 
l-fl dimensional LFD. Unless one regulates J~, the current cannot be conserved. To assure 
current conservation, it was crucial to identify the term that causes the divergence. We 
have identified this term exactly and found that it is an infinite function of the momentum 
transfer. If this infinite term is subtracted, the two LF contributions become finite as it 
must be in a conserved current. Moreover, their sum equals again the covariant result as 
expected. Still, the regularized LF contributions obtained from J~ are different from the 
ones extracted from the plus current. The differences grow with increasing binding energy. 

The covariant fermion triangle-loop (Fig. |l](a)) in 3+1 dimension is divergent if all the 
vertices are point-like. In the recent literature ||15|,|16|, the fermion- loop was regulated by 
smearing the qq bound-state vertex in a covariant manner. However, the vertex function 



used in Refs. |ll5| , p!6[| was not symmetric in the four momenta of the constituent quarks as 
we discussed in the introduction. We note that the conceptual difficulty associated with 
the asymmetry could be remedied if the cutoff prescription is used in the gauge-boson 
vertex rather than in the qq bound-state vertex. In Fig. |l|(a), we thus replace the point 
photon-vertex by a non-local (or smeared) photon- vertex S'a(A; —p)'~f'^S\{k —p'), where 
Sa{p) = A^/ (p^ — A^ + ie) and A plays the role of a momentum cut-off. Our method is gauge 
invariant, and satisfies the W-T identity |T^. Even though we have computed the covariant 



amplitude (Fig. |l](a)) with unequal constituent masses, for the clarity of presentation we 
will focus in this Section on the equal mass case, i.e. rrig = m-q = m. The basic formulas for 
the general case are given in the Appendix. 

The covariant amplitude (Fig. |T](a)) for a pseudoscalar meson is given in the equal- mass 
case by 



{p'\J^\p)=AN,g'A^ 



d^k (m^ — fc^ + p ■ p')k^ + {k'^ — rn^ — k ■ p')p^ + (k"^ — rn^ — k ■ p)p'^ 
(27r)4 D{k)D{k - p)D{k - p')Df,{k - p)Dt^{k - p') ' 

(2.2) 
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where Nc is the number of colors and g, modulo the obvious charge factor is the nor- 
malization constant fixed by the unity of the form factor at zero momentum transfer and 
the denominator factor D{k) from the quark propagator with momentum k is given by 
D{k) = k'^ — m? + ie. While the bound-state vertex is still pointlike, it satisfies correspond- 
ing relativistic bound-state equation and a binding effect of the two-body bound-state is thus 
included in our analysis. We note that one can split the denominators in Eq. (|2.2| ) into four 
terms (see below) to show the similarity between the methods of photon-vertex-smearing 
and the Pauli-Villars regularization, namely the vertex-smearing-technique shares the same 
structure of the denominators with the Pauli-Villars regularization. To investigate the issue 
of the end-point singularity 0], we present the J~ calculation in the following. 

First, the valence contribution shown in Fig. 1(b) is obtained in the range Q < k'^ < . 

By picking up the pole of the spectator quark, i.e. k~ = — — , we obtain 



{p'\J-\p) 



k+ 

,2 a4 ^ 



—AniNcg A r dk^d k±{—k^p p' +kj_-p'j_p + k± ■ p±p' — p± ■ p' j_k } 



(A2-m2)2 J (27r)4 k+{k+ - p+){k+ - p'- 

1 1 1 



E{p,A)E{p',A) E{p,A)E{p',m) E{p,m)E{p' , A) E{p,m)E{p\m) 

{2.3) 

where the energy denominator E{p,A) is defined as 

^(.■A)^^-P-- "'1.1%-/-'' . (2.4) 



As we see from Eq. ( |2.3| ), the result depends on the reference frame. However, as we will 
see later the sum of the valence and nonvalence contributions will of course be identical in 
any frame. In the present Section, we choose the frame where = ap'^ and p± = and 
the momentum fraction x is defined by /c+ = (1 — x)p'^. Then the LF valence form factor 
(see Eq. (|2.1| )) is given by 



2, 2N,g^A^ /-i f d'k^{l~x){M^ + ql)-{l + a)k^-q^ 



(2 + a) + gi/M2 7o J (27r)3 {1 - x)x^{x + 



1 



X 



1—X X I \ l~x 



1 



1—x x+a 1+a I \ 1—x x+a l+o 



(2.5) 



where M is the meson mass. The limit to the point vertex can be taken by letting A ^ cx3 
and we find that the result is finite in this frame. 

Next, the nonvalence contribution shown in Fig. 1(c) is obtained in the range p"^ < fc"*" < 

p'^ . Here, the pole is not taken as the spectator but as k^ = p'^ + " ^^l+Zp'+'' (fo'^ ^^e 

term corresponding to the Pauli-Villars particle, k~ = p'^ + ^ ^^^+1^;+'' ~''^ )- Following a 
procedure similar to the one described in the valence case, we find the LF nonvalence form 
factor to be given by 
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2N,g-'K\l + a) 



(A2 -m2)[M2(2 





a 



(27r)3a;2(l + x)(a-x) 



N{m) 



N{A) 



D0D1D2 D5D3D4: 



(2.6) 



where 



N{m) 



m' 



+ k\-{l + x) 



1 + a 



a — X 



N{A) 



1 + a 



a — X 



+ fc^ - (1 + x) 



1 + a 



a — X 



a — x 



DO 

D2 

D3 
DA 



1 + a 

M"^ + q]_ + {k± - gl) ^ m"^ + k^ 

1 + a a — X 1 + 



+ ql + (fc^ - gl) 2 , + A; , 

1 + a 



a — X 



+ 



X 



^2 _ M^ + gj ^ + (fcx - q±V ^ m^ + kj 
1 + a a — X X ' 

+ ql + (A;_L - gl)2 + fc^ 

1 + a a — X 1 + x' 

M^ + gj A^ + (fcx - gP' A^ + fcj 



- — — '-^ + 

1 + a 

2 _ M^ + gl A^ + (fcx - q±? + kj 

1+a " - + ^ • 



a — X 



X 



(2.7) 



We find from Eq. ( |2.6| ) tliat tlie two terms in tlie integrand are individually finite as long 
as the parameter A is finite. There is neither an ultra-violet divergence nor an end-point 
singularity. However, in the limit to the point vertex {i.e. A oo), we find that each term 
has not only a linear divergence in the k± integration, but also an end-point singularity in the 
X integration, which cancel each other exactly. Thus, in the point vertex limit, we find that 
the end-point singularity is completely removed even though the result is logarithmically 
divergent as it must be in the 3+1 dimensional fermion-loop with point vertices. This 
shows a striking difference from the calculation without relying on the vertex regularization 
from the beginning The critical reason for this is that the end-point singularity for the 
fermion-loop is a consequence of the bottomless nature of the Dirac sea and the vertex- 
smearing effectively provides the weighting in the Dirac sea de-emphasizing the lower part. 
The previously identified end-point singularity |^ is exactly cancelled by the identical end- 
point singularity from the term generated by the vertex-smearing corresponding to the Pauli- 
Villars particle. Therefore, in the regularized case, all the physical degrees of freedom are 
taken into account. 
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TABLE I. Our model parameters of constituent quark masses m, and cutoff values A (in units 
of GeV) used in this work. The decay constants /"y"" [p^^^] obtained by Eq. (|2]|) [Eq. ( plol) ] are 
also compared with the experimental values /^^^ [p3]. 



Meson 


rUg 


A, 


nig 


A^ 


jsym Jyasymj 


/<=^P[MeV] 


TT 


0.25 


0.90 


0.25 


0.90 


92.5 [122.2] 


92.4 ±0.25 


K 


0.25 


0.90 


0.48 


0.91 


112.5 [139.4] 


113.4±1.1 


D 


1.78 


1.79 


0.25 


0.90 


108.6 [191.5] 


< 154.9 



In addition, the sum of the valence and nonvalence contributions F(g^) = Fvai(g^) + 
Fmi.'f) is of course identical to the results obtained by other components of the current, i.e. 
either or J-*- [^,|l^.0 Also, the net result F{q^) is independent of the choice of reference 
frame. 

A similar calculation can be made for the pseudoscalar- meson decay constant /. Taking 
a non-local gauge-boson vertex, we verified again the exact cancellation of linear and loga- 
rithmic divergences in the (two-point) fermion loop. The result for the equal mass case such 
as the pion is given by 

N^gmk^ /•! , , rC(m,m)C(A,A)' 
/ dx log 
/o 



/ 



4v^7r2(A2 
where C {mi, 1712) is given by 

C(mi,m2) = 



m 



2\2 



C(m,A)C(A,m) 



(2.8) 



x(l - x)M^ 



'1 



XTTln 



(2.9) 



If the asymmetric qq bound-state vertex-smearing is used as in Refs. |]T5| , p^ , the above result, 

NcQmA'^ , ^ C{m,m) 



Eq.( |2.8|) , is replaced by 
/ = 



dx log 



C(m,A) 



(2.10) 



4v^7r2(A2 -m2 

We note that the logarithmic divergence is not completely cancelled in the unequal mass 
case if the minus component of the weak current is used with the asymmetric qq bound-state 
vertex-smearing. The non-local gauge-boson vertex used in this work does not suffer from 
such an incomplete cancellation of divergences no matter which component of the current 
is used. 

In the next section, we present the numerical taxonomy of the vr, K, and D meson 
electromagnetic LF form factors choosing various reference-frames, as well as the values of 
the decay constants f^,, fx and fjj. 



III. NUMERICAL RESULTS 

In Table |I|, we present our model parameters such as the constituent quark masses (m), 
and the cutoff values (A), as well as the decay constants that we calculated here and com- 
pared to the experimental data [E4|. The meson masses (M) are taken as the experimental 



2ln the electromagnetic form factor, there is no zero-mode contribution in J± current. 
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TABLE II. The kinematics in the reference frames used in this work, where k = Q^/2M and 



h = (cos (j), siiKp). 



Kinenicitics 


Target rest frame 


Breit frame 


u 1 vv irame 




K + Qy^l + k/2Mcos6' 


+Q cos 6 







K - Qy/1 + k/2Mcos6' 


—Q cos 9 




q± 


Q^/l + k/2M sin 6*71 


Q sin 9h 


Qn 


p+ 


M 


v/M^ + QV4 - g+/2 




p- 


M 


VM^ + QV4-g-/2 




p± 





-9l/2 






values Our model parameters have been chosen to fit both the charge radii and the 

decay constants to the experimental data well, although the available data for the decay con- 
stant of the D-meson is only an upper limit. We also compare in Table | the decay constants 
from our symmetric non-local gauge-boson vertex (/^y™) with those from the asymmetric 
qq bound-state vertex (/'^^y™) smearing case. However, we note that our calculation here is 
limited in value because the zero-range approximation is used for the bound-state vertices. 

Three different reference frames were considered: the Drell-Yan-West (DYW) frame, the 
target-rest frame (TRF) and the Breit frame. The corresponding kinematics is given in Ta- 
ble The DYW frame has gained some popularity in deep-inelastic scattering calculations 
because in that frame g"*" = identically. This frame can be obtained by taking the a ^ 
limit in the frame presented in section II. In the limit the frame presented in section 

II coincides with the target-rest frame with ^ = 0. 

Now we comment on our results shown in the figures below. (In all these figures, we use 
thick solid line for the covariant form factor, thick dashed line for F^^, thick dot-dashed line 
for F^, thin dashed and dot-dashed lines for the corresponding minus LF form factors.) We 
show in Fig. the results of our numerical calculations using the DYW frame and compare 



with the experimental data p8|-p3| Our total results, represented by the solid line in each 
figure, are also in very good agreement with the experimental data of the pion and kaon 
form factors, respectively. 

In the DYW frame F^^ vanishes identically. Remarkably, we find that the nonvalence part 
of the minus current, which in this reference frame coincides with the zero- mode contribution, 
makes a very important contribution to the total form factor and may even dominate over the 
valence part in the whole Q^-range considered. There are quantitative differences between 
the results obtained for the different mesons, which are to large extent due to the difference 
in binding: the tighter the binding, the more important F^ becomes. We checked the 
binding effect in the case of the pion. By varying the quark mass alone from the realistic 
value of 0.25 GeV to 0.07007 GeV, thus lowering the binding energy to 0.1% of the pion 
mass, we found that the value of at = was increased from 6% to 69%. Still, for 

> 0.3 GeV2 

-^vai is lower than F^^. This indicates that for the larger values of Q^, the 
relativistic effects can still be large. 

The frame dependence of the different components of the current can be studied by 
comparing the results of calculations in different frames and at different values of the polar 
angle 6. It is worth mentioning that the results must be independent of the azimuthal angle 
0, because rotations about the 2;-axis are kinematical transformations. (We used this fact 
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check on the correctness of our codes.) 

In Fig. we show the resuhs of our numerical calculations in the Breit frame at 6' = 0, 
which differs from the target-rest frame at ^ = by a boost in the z-direction only, so the 
results are identical in these two frames. The first thing we notice is the great difference 
with the DYW frame. Now there is a sizeable contribution from F^, which dominates at 
higher values of the momentum transfer: at = 0, coincides with the covariant form 
factor. It crosses at some value of Q^, the crossover point being smaller for larger binding 
energy. It is of special interest to separate the instantaneous part, i.e. the contribution to 
F~^ from diagrams with one internal instantaneous propagator. They are given in the figures 
by crosses (x). It turns out that they give a very large contributions to the pion form factor 
but become negligible for the heavy Z)-meson case. 

Turning to the minus current we see that in this reference frame the relative importance 
of F^ becomes more prominent than in the DYW frame for the more tightly bound mesons 
TT and K. In the case of the pion the dominance of F^ is so strong, that it can hardly be 
distinguished from the covariant form factor. The figures (in Fig. |^) suggest that at = 
the values obtained are frame independent. This is indeed the case, as we found in our 
calculations and can be understood as the equality of the form factors in the long-wave- 
length limit. 

A first glance at the angle dependence is given in Fig. ^. In both the target-rest frame 
and the Breit frame for 9 = 7i/2 the transferred momentum is purely transverse. As the 
pure Lorentz boosts in the transverse direction, Ki and K2, are dynamical in LFD, we must 
expect a strong frame dependence and this is indeed what we find. As g"*" = in the Breit 
frame for 6 = 7r/2, the part F^ vanishes identically. Still, the DYW frame results differ 
for the minus current, which is clear from a comparison with Fig. 0. The reason is that 
they are not connected by a kinematical transformation. For the same reason the zero-mode 
contributions differ in the two frames. In the DYW frame it is very close to the total form 
factor, while in the Breit frame it even overshoots the covariant form factor at large values 
of by a factor of almost 2. In this connection we want to mention the work of Frederico et 
al. [l^, who performed a calculation similar in spirit to ours, but claimed that in the Breit 
frame at 6 = 7r/2 is always very close to the covariant form factor. This discrepancy 
was due to the difference of their definition in the Breit frame at = 7r/2, where they 
removed the term that is odd under the transformation p~ — k~ — ^ —{p^ ~ ^ ) their 
definition ||2^. However, our definition of F~^i is general for any angle 9 and our results on 
^-dependence are smooth as shown in Fig. ^ 

We show the systematics of the angle dependence for the case of the pion only. Our 
results are depicted in Figs. ( The thick solid line for 6* = 0, thin dotted line for 

9 = 7r/4, the thick dashed line for 9 = 7r/4, the thin dot-dashed for 9 = 3tt/4, and the 
the thick long-dashed line for 9 = n.) One sees immediately that the angle dependence is 
smooth but can be very strong, both for the valence and the nonvalence parts, calculated 
from either the plus or the minus current. 

One might try to exploit the angle dependence to optimize the calculation of the form 
factor in a noncovariant framework. However, as the figures above clearly show, there is no 
value for the angle 9 where both F^ and F^^ are negligible, or even suppressed, compared 
to the valence parts, for all values of Q^. On the contrary, as the values of the form factor 
components at = are frame independent, we can be sure that F~^ must be very 
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important for an important part of the range. 

As a summary of what we found concerning the angle dependence we show in Fig. ^ 
the complete angle dependence for two values of the momentum transfer, = and 0.1 
GeV^. All curves in this figure are clearly smooth and demonstrate the fact that there is no 
preferred value for the polar angle Q. 

IV. CONCLUSION AND DISCUSSION 

In this paper, we have analyzed all the components of the current quantized on the 
light-front to compute the electromagnetic form factors of pseudoscalar mesons with spin- 
1/2 consituents. Since our aim in this work was to analyze the taxonomy of the triangle 
diagram, we did not choose any particular electromagnetic gauge but just presented the 
equivalence of the physical form factor F{q^) in any choice of the current component. The 
divergence appearing in the 3+1 dimensional fermion-loop calculations was regulated by the 
covariant vertex-smearing-technique. Performing the light-front calculation, we verified that 
the vertex-smearing-technique is similar to the Pauli-Villars regularization. In the J~ com- 
putation, we find that the critical smearing effect persists even in the limit to the point vertex 
because the end-point singularity existing otherwise is completely removed once the limit is 
taken at the end of the calculation. If the limit is taken at the beginning of the calculation, 
however, we have alreay shown that the end-point singularity in the nonvalence contribution 
leads to an infinitely different result from that obtained by the covariant Feynman calcula- 
tion . Our taxonomical analysis demonstrated that each individual contribution, whether 
valence or nonvalence, is finite regardless of which component of the current is considered. 
However, we stress that each contribution depends on the reference-frame even though the 
sum does not. Of course, the invariance of the sum ensures the current conservation. Also, 
the zero-mode contribution should be distinguished from the instantaneous contribution as 
we have numerically estimated the differences. 

lYiom. our numerical calculations we can conclude that to obtain agreement with the 
covariant form factor one needs both the valence and the nonvalence parts. For tightly 
bound states the nonvalence parts dominate in an important part of the range of Q^-values 
that we studied. It is natural that this result runs counter to nonrelativistic intuition, 
which says that the valence parts should dominate, because the tightly bound states are not 
expected to be non-relativistic. In Ref. it was demonstrated in the 1+1-dimensional case, 
that by weakening the binding and F~^^ approximate the covariant form factor more and 
more closely. Here we found that indeed for less tightly bound state the valence parts come 
closer to the covariant result, but even at = gives only 69% of the covariant form 
factor. 

If two reference frames can be connected by a kinematical Lorentz transformation the LF 
form factors calculated in these two frames must be the same. Otherwise they must differ. 
We found that the angle dependence within a chosen reference frame is always smooth, 
although it may be very strong. For some values of the polar angle 6 two reference frames 
may be connected by a kinematical Lorentz transformation, e.g., the target-rest frame and 
the Breit frame for ^ = 0. On the other hand, the DYW frame and the Breit frame for 
6' = are not connected by a kinematical transformation and the form factor components 
consequently do not coincide. More details of the frame-dependence can be found in Ref. ||26| . 
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While the calculations carried out in this work is more realistic than the case of 1+1 dim 
with the point vertex, they are still semi-realistic as the model uses a bound-state rather than 
a confined-state. Moreover, we used for the vertices those obtained from a Bethe-Salpeter 
equation using a contact interaction which are of zero range and may emphasize the 
importance of the nonvalence contributions rather differently from a more realistic model. 
Thus, there is still much room for extending our model towards a more realistic model. 
However, the essential conclusions about the similarity to the Pauli-Villars regularization 
and the cancellations of both ultra-violet divergence and the end-point singularity remain 
intact. Nevertheless, the numerical results may differ from a more realistic model calculation. 
This point is presently under investigation. 
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APPENDIX A: UNEQUAL-MASS CASE 



If we number the momenta of the internal hnes of the fermion triangle as ki (spectator), 
^2 = P + ki, and k^ = p' + ki (struck quark) respectively, and the corresponding masses 
as mi, m2, and then we find for the trace appearing in the numerator of the covariant 
integral, the expression 



-(m2m3 - k2± ■ k3±)k'^ + {rriims - k2± ■ k3±)kl^ + (mim2 - k2± ■ k3i_)k^ . (Al) 



The valence diagram is obtained if one calculates the integral over ki by closing the contour 
around the pole corresponding to putting the spectator ki on shell. This corresponds to the 
following values for the minus-components of the momenta to be used in the expression for 

ki = k^ =p- + fcf, fcg- = p'- + /cf. (A2) 

To obtain the nonvalence diagram one closes the contour around the pole corresponding to 
k^. Then one gets for the part of the propagator with the mass rris 

K= '\"t'"' > K = -p'' + K^ K =p' -p'' + K- (A3) 

For the part of the propagator with the cutoff A' the same formula can be used, but with 
m3 replaced by A'. 

As the smearing we use affects only the denominator parts of the propagator of the struck 
quark, the replacement of by A' in T'^ occurs only in the minus components. 
The energy denominators are easy to find. One obtains for the valence parts 

D{k - p') = p" + ^L±i!i _ Mi±z!i. (A4) 

kl % 

The nonvalence part has the same denominator D{k — p'), but D{k — p) is changed to 



k^ L ~^ ^3 "I" 



D'{k -p)=p--p'-+ ' - '\^ ' . (A5) 



The energy denominators D\ and D'^, are obtained by the substitutions 1712 A and 
m3 — > A', respectively. 

In the case where mi = 1712 = and A = A' the final formulas are much simplified. 
The explicit formulas in the main text are valid for the equal-mass case. The general case 
can be easily constructed from the expressions ( [Al| - |7^ ). 

In the explicit formulas we use the notation 

i,t^ i-f^ ^ 1 ^ k+ -p^ p'+ -p+ 

k^ = —kl, X = — (valence), x = ; (nonvalence), a = ; . (Ao) 

p+ p+ p+ 
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APPENDIX B: ANALYTICITY 



/^From the covariant expression for the amphtude one can prove that the form factors are 
analytic functions of . This proof is not vahd for the LF-time-ordered amphtudes. One 
can, however, expand the expressions that one obtains for the different parts of the form 
factor in terms of Q and determine the Q-dependence at small values of . It turns out 
that the results depend on the kinematics: it matters in which reference frame one does the 
calculations. 

We use the fromulas for the momenta given in Table || and expand the trace and the 
energy denominators in powers of Q. From this expansion the analyticity properties of the 
amplitudes follow. We have numerically verified that a blow-up of Figs. for small 
illustrates the analyticity properties of the amplitudes discussed in this Appendix. The three 
different reference frames are discussed consecutively. 

1. Target-rest frame 

We discuss the valence part of the plus-current in detail; the other LF form factors can 
be treated in a similar way. First we expand the momenta. 

g± = gcos^ + C(g2), gl = Qsin^n + 0(Q^), 

p'± = M + Qcos^ + C(Q2), p1 = Qsin^n + C(g^). (Bl) 

The trace consists of a piece that is independent of Q and a piece that consists of two 
parts, one that is proportional to Qki± ■ -hsinO and a part proportional to QcosO. The 
denominator D{k — p') has a similar behaviour, but D{k — p) is independent of Q. 

Upon integration over ki^ the terms proportional to • nsin6' vanish. Consequently, 
the valence part has the small-Q behaviour 

^vti~i^vli° + i^vti'Qcos^. (B2) 

The other cases, nonvalence plus-current, valence and nonvalence minus- current, show 
the same pattern. There are pieces independent of Q, parts with the Q-dependence Qki± ■ 
nsinO and ones proportional to Qcos6. As the pieces proportional to ki±- n vanish upon 
integration, all the components of the form factor show a behaviour similar to Eq. ( |B2[ ). 
So, only for 6^ = 7r/2 do we find that the LF time-ordered amplitudes calculated in the 
target-rest frame are analytic in Q^. 

2. Breit frame 

In the case of the Breit frame we can follow the same line as in the case of the target-rest 
frame. The only difference is that now no terms of the form Qki± ■ nsinO appear. As those 
terms give no contribution to first order in Q anyway, this does not alter the result: also in 
the Breit frame the components of the form factor have the same small-Q behaviour as in 
Eq. (il). 
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3. Drell-Yan-West frame 



In the Drell-Yan-West frame there is no angle dependence. However, there is a term linear 
in Q. It is proportional to Qn ■ ki^. Of course, this term also vanishes upon integration over 
ki^. Therefore, in the Drell-Yan-West frame we find no term linearly dependent on Q, so 
the amplitudes are analytic in Q^. 

4. Zero mode 

The zero mode is defined as the contribution to the nonvalence amplitude that survives 
the limit q'^ — >■ 0. It is easy to see that only the minus-current can have a zero-mode part, 
because the integral defining has an integrand that remains finite when goes to 0. 
However, the integrand defining diverges when the limit g+ — is taken. In order to 
determine the limit, one may expand the integrand in powers of q^. As q~^ is taken to be 
zero in the Drell-Yan-West frame, the analysis cannot be done in that frame, but it can 
most easily be carried out in the Breit frame. 

The algebra being straightforward but tedious, we shall not give the details. Rather 
we quote the final result. If we consider the smeared case, the most divergent part of the 
integrand for F^ has the behaviour As the integration over kf ranges from —p'^ to 

— p"*", we can scale the integration variable as in Eq. ([A6| ) and obtain an integral over x from 
to 1. The Jacobian being g"*", it cancels one factor g+ in the denominator, so we only need 
to show that the leading term vanishes to prove that the zero-mode contribution is finite. 

Upon carrying out the algebra we find the following result. The leading part is propor- 
tional to 

x\kl^ + ml){kl^ + A'2) - (1 - xf{kl^ + ml){kl^ + A^). (B3) 

This part of the integrand does not vanish, because A;2± = P± — k± ^ k-^^ = f>'± ~' k±. 
However, if we take A' = A and m2 = m^, being the mass of the struck quark, and as 
k2±, ^3±, and k± differ only by a constant vector, this function vanishes after integration 
over X from to 1 and k± over the whole of R2. So we see that the coefficient of the 
contribution proportional to l/q^ vanishes. The remaining part, the piece that survives the 
limit g"*" — s> 0, is the zero mode contribution. In the Breit frame as well as the target-rest 
frame it is well-defined and finite. 
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FIG. 2. Pion, K, and D meson form factors in the Drell-Yan-West (DYW) frame compared 
with experimental data for pion [p!^22| and kaon |23]. The part vanishes identically in this 
reference frame. 
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FIG. 3. Pion, K, and D meson form factors in the target-rest frame (TRF) or the Breit frame, 
both at 6* = 0. 
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FIG. 4. Pion form factor in the target-rest frame and the Breit frame at 9 = 7r/2. 




FIG. 5. Pion LF form factors and in the target-rest frame for five different values of 
the polar angle 6. 
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FIG. 6. Pion LF form factors and in the Breit frame for five different values of the 
polar angle 9. 




FIG. 7. Pion LF form factors F^^^ and F^^ in the target-rest frame for five different values of 
the polar angle 9. 
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FIG. 8. Pion LF form factors F^^^ and F^v the Breit frame for five different values of the 
polar angle 9. 
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FIG. 9. Systematics of the angle dependence of the LF form factors at = 0.1 GeV^(thick 
lines) and GeV^(thin lines) in the Breit frame. 
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